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| Abstract. Precise mechanisms by which Active Galactic Nuclei (AGN) re- 

■ ceive their gaseous fuel is still a mystery. Here I draw attention to the extra 

' ordinary star formation event that took place in the central ~ 0.5 parsec of 

our Galaxy. The most reliable explanation of the event seems to be that two 
somewhat massive nearly co-eval gaseous disks failed to accrete on Sgr A* , the 
super-massive black hole (SMBH) in our Galaxy, and instead cooled down and 
gravitationally collapsed, forming the stars observed now. This emphasises that 
star formation must be an important part of AGN feeding puzzle. I also discuss 
a model in which stellar winds create the observed obscuration of AGN. These 
winds are cold, clumpy and dusty, as required by the observations, but they are 
Compton-thin unless wind outflow rate is highly super-Eddington. This argu- 
ly-j | ment is in fact a general one, independent of the wind driving mechanism. I 

thus suggest that winds may be important for optically thin absorbers, and that 
a better model for optically thick AGN obscuration is a warped accretion/star 
forming disk. 
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1. Introduction 
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In X-ray binary systems, we know quite well the origin and properties of gas 
arriving at the outer edge of an accretion disk. In contrast to this, we know noth- 
ing for a fact about how gas trickles down into the inner parsecs of AGN. This 
is because observations of these regions are tremendously difficult: one needs to 
resolve very small angular scales faced at the same time with the contaminating 
emission from few tens of gravitational radii of SMBH. Despite all this, there is 
now a significant progress in such observations ( Davies et al. 20061 ) . 



The Galactic Center (GC) is our most suitable laboratory for these stud- 
ies due to its proximity and the fact that Sgr A* is in fact a very dim source, 
producing bolometrically only ~ hundreds of Solar luminosities. As a result, 
young massive "He-I" stars dominate the power output of the central parsec 
of our Galaxy. Many of these stars are now resolved to be located in a very 
well defined and rather thin stellar disc that rotates clockwise as seen o n the 
sky (|Levin & Beloborodov 19931 : iGenzel et al. 20031 : iPaumard et al. 20061 ) . The 



rest of the stars can be arguably classi fied as a second mo r e diffuse disc or a fea- 
ture that rotates counter clock-wise (IGenzel et al. 20031 : IPaumard et~ai~ 2006: 



Lu et al. 20061 ). 



Below I shall argue that the best interpretation of data on these young 
massive stars is star formation inside an accretion disk(s) that existed in the 
inner parsec just before the stars were born. I shall discuss some observational 
and theoretical issues pertaining to formation and mass spectrum of these stars. 
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I will also present results of numerical simulations of star formation in a massive 
gaseous disc around Sgr A*. 

Using Sgr A* observations as a launch pad for a discussion on the state 
of gas arriving in the inner parsecs of AGN, I will then argue that in general 
AGN disks are star-forming as well. I will then present numerical simulations 
of stellar winds escaping from a stellar disc left by a star formation event. Such 
winds might provide some obscuration of AGN, as required by observations. 
The results of the simulations however show that such outflows are unable 
to provide Compton-thick obscuration unless winds carry significantly super- 
Eddington outflow rates. A simple analytical estimate shows that this result 
is general for any wind-driving mechanism. I thus argue that, while certainly 
present and important for the working and appearance of AGN, outflows are 
not responsible for Compton-thick obscuration. Instead, I propose that the role 
of the Compton-thick absorber is played by a geometrically thin but strongly 
warped accretion/star- formation disk. As an example, gravitational torques due 
to a non-spherical stellar potential are demonstrated to be effective in producing 
strongly warped gaseous discs. 



2. Observational constraints on the young stars in Sgr A* neighbour- 
hood 

2.1. The in-situ origin of the stars 

There is nearly a hu ndred early type stars in the central parsec of our GC 
( Paumard et al. 20 06). Finding these young stars there, as close as ~ 0.03 



pc away from a SMBH, was a real surprise. "Standard" models of star for- 
mation are not easily applicable here due to a huge tidal field of the central 
object at R = 0.1 pc distances from Sgr A*. The required gas density is 
nu > 10 n cm~ 3 (i?/0.1pc)- 3 , i.e., exceedingly large compared with gas densi- 
ties in galactic molecular clouds. 

It thus has be en suggested th at the stars may have been formed at a distance 
of tens of parsecs ( Gerhard 200 ll ). avoiding the need for the excessive gas density 



prior to star formation, in a massive star cluster. The cluster's orbit would 
then decay through d ynamical friction with the background stars. However, 
( Paumard et al. 20061 ) find a well defined outer edge of about ~ 0.5 parsec to 



the observed distribution of young massive stars in the GC. This is difficult to 
understand if stars migrated into the inner parsec from outside, as the disrupted 
cluster should leave a trail of stars behind it. 

If the cluster were very strongly mass-segregated, i.e., that massive stars 
were present only in its core, then the trail could be composed of low mass 
stars only. These stars could escape detection in the near infra-red. This model 
seems to be challen ged by X-ray data. R ecent calibration of X-ray properties 
of T Tauri stars by ( Preibisch et al. 20051 ) revealed that they are very bright in 



X-rays, i.e. some three orders of magnitude brighter than they are on the main 
sequence. They would be observable even at high extinction regions such as the 
GC. These young stars also display giant X-ray flares, with some becoming as 
bright as ~ 10 33 erg s _1 . ( Navakshin Sunvaev 20051 ) used these new X-ray 



results for young stars, and observations of the GC reported in (? ), to show 
that the inner ~ 10 parsec could not be hiding more than 10 4 or so young low 
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mass stars. This is insufficient: some 10 5 — 10 6 Solar masses of young stellar 
mass is needed to make the cluster heavy enough for it to sink in during the 
short lifetime of the young massive stars. 

2.2. Top-heavy mass spectrum of the in-situ star formation event 

Because Sgr A* is so dim in X-rays, it turned out possible to push the X- 
ray constraints further. Baganoff et al. (2003) described the properties of the 
unresolved X-ray emission near Sgr A* in detail. Nayakshin & Sunyaev (2005) 
used these results to deduce that the area most densely populated by massive 
young stars, i.e., the inner R < 0.2 parsec could contain no more than ~ 10 3 
Solar masses of low-mass YSO (M < 3 Solar mass). This is interesting since a 
factor of 10 o r so more would be ex pected if the IMF were the "normal" galactic 
one, such as (|Miller fc Scab 19791 1. 



2.3. Total stellar mass created 

An important question to ask is the total mass budget of the star formation 
event. The best age estimate for the young stars in the G C is ~ 6±1 million year s 



and their combined mass is a few thousand Solar mass (jPaumard et al. 2 006). 
Due to the age of this stellar population, stars more massive than ~ 50 — 60 Solar 
masses seem to be absent. Due to the top-heavy nature of the mass spectrum of 
the young stars discussed above, it is then possible that the initial star formation 
event created a star cluster far more massive than the one we observe now. 



( Navakshin 2005 ; Navakshin Sz Cuadra 20051 ) proposed a way to constrain 



the total stellar mass using observational constraints on stellar orbits. Suppose 
that the two stellar systems were created infinitely thin and flat, as they would 
be in the simplest self-gravitating disc scenario (see §3 below). With time, an 
isolated stellar disc will thicken due to internal A-body heating, and two stellar 
discs will warp each other due to their non-spherical gravitational potentials. 
Both of these effects are stronger the more massive the stellar discs are. Now, if 
the initial systems wer e not thin and flat, th en the thickening and warping will 
be even quicker (e.g., Navakshin et al. 20 06). Therefore, demanding that the 



model orbital configuration fits two stellar planes no worse than the real stars 
do at the present moment, we can arrive at a constraint on the total mass of 
these discs. This yields the limit of around 10 4 solar masses for the total masses 
of each of the stellar discs. 



Theoretical models 



It has been known for a long time that massive accretion discs can be grav- 
itationally unstable and form stars if they are able to cool efficiently (e.g., 
Paczvnski 19781 ; Gammie 2001 ; Rice et al. 20051 ). In particular, the disc has 
to be both massive enough and be able to cool fast enough: Toomre parame- 
ter Q < 1, and i coo i < 3S7 1 , where O, = GMbu / R 6 ■ Vertically integrated 
marginally star-forming disc models have been consid ered recently with applica- 
tion to Sgr A* in mind ( Levin 20061 ; Navakshin 20 06). Pleasingly, these models 
have a number of features consistent with the observations. The inner edge of 
the self-gravitating region, i.e. the closest to Sgr A* the stars can be born in 
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this scenario, is about 0.03 pc, consistent with the distribution of the disk-stars 
( Paumard et al. 20 06). Further, the models predict that all the gaseous disk 
mass would go into building up stars. The stellar mass is then estimated to be 
at ~ 10 4 Solar masses, which is not inconsistent with the observations. 

The issue of the top-heavy mass function is more interesting. The estimates 
for the fragmentation mass, Mf rag , i.e. the mass of the first bound fragments, 
are certainly below 0.1 Solar mass. Hence, if disk were to promptly collapse into 
clumps of mass of this order, one would expect low-mass stars to dominate the 
mass spectrum of collapsed objects. However, in reality, collapse of gas clumps 
may not be dynamical (as assumed by the simple model), but gradual, regulated 
by cooling and clump rotation. Further, the clumps will have sizes comparable 
to the dimensions of the "first cores", i.e. -Rr.iump ~ 5 AU. The clumps are 
actually likely to merge and grow by agglomeration ( Levin 20061 ). Therefore one 
expects that Mf rag is a strong under-estimate of the actual stellar mass. Further, 
Nayakshin (2006) noted that low-mass proto-stars born with M ~ Mf rag may 
grow very massive by gas accretion if disc fragmentation is slow. In particular, 
stellar accretion luminosity, produced by newly born proto-stars, turns out to 
be sufficient to heat the disc up, increasing the Toomre Q-parameter of the 
disc above unity and hence making the disc stable to further fragmentation. 
Numerical simulations (see below) do bear this prediction out. 



4. Numerical models of star-forming discs 

We use the SPH/N-body code GADGET-2 (ISprineel 20051 ) to simulate the dy- 
namics of stars and gas in the (Newtonian) gravitational field of Sgr A*. Details 
of our method will be reported in Nayakshin, Cuadra & Springel 2007. Radia- 
tive cooling of the disc is treated with a simple locally constant cooling time 
prescription, i coo i = /3/0, where (3 is a constant of order unity. 



Circular gas discs without feedback. Left panel of Figure Q] shows a snap- 
shot of a run with (3 = 3 well into the non-linear stage, at time t ~ 10 4 years, 
when more than half of the gas was already turned into stars. The initial con- 
dition is a gas disc of mass 2 x 10 4 Solar mass in Keplerian circular rotation 
around Sgr A*, extending from 1" to 4" (1"~ 0.04 pc). The right panel shows 
the same simulation in a slightly earlier time but now zoomed into a smaller 
patch of the disk to emphasise the disk structure. 

Star formation is fastest in the innermost region, where most of the gas is 
already depleted. At the end of the simulation essentially all the gas is turned 
into stars. Unlike "normal" star formation, feedback from massive stars will not 
be very effective in blowing the gas away via radiation pressure or winds as the 
stars are within the deep potential well due to Sgr A*. The stars thus steal the 
majority of SMBH's dinner. 

Interestingly, gravitational heating generated by stars (scattering off each 
other and interacting with gas) is sufficient to heat the disc up above its pre- 
star-formation value, slowing down and even shutting off fragmentation at later 
times. The effect is more pronounced the longer the cooling time. Figure [21 left 
panel, shows fragmentation rate in Solar masses per year as a function of time 
for simulations that differ only by the cooling time parameter /3, as labelled in 



5 




the caption. The longer the cooling time, the sooner disc fragmentation slows 
down and eventually stops as the disc becomes hotter. The total mass budget of 
the simulation is fixed by the initial condition. Quite logically then, the longer 
the cooling time the more top-heavy the mass function of the stars born in 
the simulation becomes. The IMF of the stars from these three simulations is 
presented in Figure [21 right panel. 



The role of feedback on the fate of the discs and the IMF. As already 
mentioned, Nayakshin (2006) suggested that the accretion luminosity of young 
stars within a star-forming disk might be strong enough to warm it up and 
significantly slow down its fragmentation. We have made several numerical runs 
with stellar accretion feedback included. As expected, we found a very strong 
reduction in the fragmentation rate, and an increase in the disk lifetime defined 
as the time it takes for a half of its initial mass to be turned into stars. We 
could not reach a conclusive statement on the influence of the stellar accretion 
luminosity feedback onto the stellar mass function as the simulations could not 
be run sufficiently long. 

The fact that disk lifetime increases strongly in the simulations is encour- 
aging in terms of satisfying SMBH growth needs. If these results could be 
extrapolated into a much heavier accretion disc case, the time scale for angu- 
lar momentum transfer mediated by self-gravity could become shorter than the 
disc lifetime. Unfortunately, this appears unlikely. In the simulations we used 
the simple (3 cooling time prescription. Quite a general analytical argument 
shows that, in order to feed AGN at rates of order Eddington accretion rates, 
the required levels of energy feedback from star formation (and in fact of any 
extra disc heating at all) a re so large that t hey strongly contradicts the observed 
spectra of a typical AGN (jGoodman 20031 ). 
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Figure 2. Left panel: Disc fragmentation rate in the tests with (3 = 0.3 
(thick solid), (3 — 2 (dashed), and j3 = 3 (thin solid). The shorter the cooling 
time, the faster disc fragmentation proceeds. Right panel: Mass function of 
stars ("IMF") formed in the simulations with /3 = 0.3 (solid), (3 = 2 (dotted), 
and (3 = 3 (dashed). The dot-dashed power-law is the Salpeter IMF. 



Star formation in eccentric discs The second of the observed stellar discs in 
the inner parsec of our Galaxy is a not so well defined diff use feature with stella r 
orbits that appear to be rather eccentric (e ~ 0.5 — 0.8 Paumard et al. 2 006). 
Since N-body relaxation time scale is about a Gyr, these stars must have been 
born on already eccentric orbits. We performed a series of tests aiming at 
establishing whether star formation can take place in eccentric discs. The results 
are rather similar to the runs with circular gas discs. 

We think that the best case scenario for the GC stellar populations is one 
where the clock-wise disc stars formed out of a more massive circular gaseous 
disc, whereas the counter clockwise stars were born out of an eccentric stream 
of gas. Orbital precession could then turn an initially flat eccentric disc into a 
geometrically thicker one by now. 



5. Star- forming winds and AGN obscuration 



The talk by Moshe Elitzur (jElitzur 2 006) summarises the observational need for 



a dusty AGN absorber and reviews important theoretical issues in that topic. It 
is proposed that a dusty clumpy outflow might be the best theoretical explana- 
tion for the absorber. Here I discuss related ideas. In particular, star formation 
on sub-parsec scales, discussed in previous sections of this paper, will naturally 
lead to stellar outflows. Provided the mass outflow rate is large enough, winds 
will be dense and thus cold, and quite likely dust-rich. In this case one might 
expect these outflows to play an important role in the obscuration of AGN. 

SPH numerical modelling of stellar wind accretion on Sgr A* was recently 
developed by us (Cuadra et al. 2005,2006). Nayakshin & Cuadra (2007) used 
that numerical approach to model stellar wind obscuration, concentrating on 
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higher mass outflow rates than in Sgr A* presently, as that might be expected 
in a bright AGN environment. 

In the simulation presented here, stellar mass loss rates are 2.5 x 1CP 4 
Solar mass year -1 per star. In total, we have 200 mass shredding stars, thus 
amounting to the mass loss rate of 0.1 Solar mass per year, which is a factor of 
few super-Eddington for Sgr A* mass of M = 3.5 x 10 6 Solar mass. The stars 
are situated in a flat circularly rotatingKeplerian disk of geometrical thickness 
H(R) = 0.1R, where R is the radiuqj. The disc inner and outer radii are 
-Rin = 1.5 and R OVLt = 8, respectively. 

A snapshot of the simulation is shown in Figure [3] at time t ~ 2200 years 
after the beginning of the simulation. While there is no true steady state for a 
system of a finite number of moving stars, the snapshot is fairly typical of the 
morphology of the stellar wind. The face-on view (left panel) shows that some of 
the shocked wind managed to cool down and formed a small-scale disc. Smaller 
escape velocity and less frequent shocks at the larger radii allow direct escape 
of both the fast diffuse and the slower cooler clumpy winds. Due to the final 
extent of the stellar disc and a projection effect, the wind morphology reminds 
an "X" -shape (Fig. [31 right panel). 
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Figure 3. Face-on (left panel) and edge-on (right panel) views of the simu- 
lation domain for the flat stellar disc configuration. Black asterisks show the 
location of the wind-producing stars used in the simulation. From Nayakshin 
& Cuadra (2007) 

The focus of our attention here is on the obscuration properties of these 
winds, and not on that of gaseous discs. The left panel of Figure H] shows the 
obscuring column depth of the winds as seen from the SMBH, excluding the gas 
with radial distances R < 1.6 from Sgr A*, which removes most of the gaseous 
disc. Notice the very irregular patchy structure of the (brighter) optically thicker 



J The unit of length used here is 1 arcsecond, which corresponds to about 1.2 x 10 17 cm or 0.04 
pc at the Galactic Centre distance of 8 kpc. 
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regions. The contrast between those and neighbouring less dense patches of sky 
is frequently a factor of 10 or more. The dotted pattern at the cos 9 = plane 
are the dense regions of stellar winds immediately next to the stars. 
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Figure 4. Column depth through the wind (left) and isotropic mass outflow 
rate (right; in units of Solar mass year -1 ) for the simulation shown in Fig. [3] 
as seen from the SMBH. Note the large variations of the obscuring column 
depth over small angular scales. The gaseous inner disc seen in the left panel 
of Fig. [3] has been excluded. 

The right panel of Figure H] shows the isotropic mass loss rate along the line 
of sight, defined as 

M= Jd ^ 4nr2pVR (1) 

where the integral is taken along the line of sight defined by a given 6 and 4>. 

Now, while the obscuration pattern in the left panel of Figure U] may cer- 
tainly be relevant to the observational need for a cold dusty absorber, there 
appears to be a serious limitation for this model. The average column depth 
is not large enough, i.e. it is Compton-thin. Ho wever, observations require a 
large fraction of AGN to be Compton-thick (e.g., ISazonov & Revnivtsev 20041 : 
iGuainazzi et al. 20051 ). To satisfy this requirement, the mass outflow rate, al- 
ready super-Eddington by a factor of ~ 3, needs to be pushed another factor of 
> 5 higher. 

This is an entirely general point. Considering clouds with average column 
depth S c outflowing at speeds about t he escape velocity, and num bering iV a per 
line of sight on average, one arrives ( Navakshin fe Cuadra 20071 ) at the mass 
outflow rate of 

Mq / r. \V2 

Mwind ~ 15- 



n 



M 



1/2 



, . (2) 
year \ 1 pc/ 

where Mg is the SMBH mass in units of 1O 8 M0, and r t is the "torus size". 
Now, the Eddington accretion rate is AirGMrrip/ecaT ~ 2MQyear _1 Ms. The 
outflow is Compton-thick when iV a E c > 1. Thus the required wind mass loss 
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rate (equation [2D is an order of magnitude higher, typically, than the Eddington 
accretion rate for a Mg ~ 1 object and a parsec-scale torus. 

Considering a specific case of the local obscured AGNs studied by Guainazzi et. 
al. (2005), we note that the bolometric luminosities of these objects in the in- 
frared, X-ray and optical bands are in the range L ~ 10 43 — few x 10 44 erg/sec, 
which implies SMBH accretion rates of "only" ~ O.OIM0 year -1 for the stan- 
dard radiative efficiency. Hence if the obscuration of the optically thick objects 
in that sample were provided by the winds, we would conclude that the SMBH 
accretion process must be very wasteful, with ~ 100 — 10, 000 times more mass 
flowing out of the inner parsec than accreting on the SMBH. It would also re- 
quire a very high mass influx into the inner parsec to sustain such winds. Given 
the difficulty of delivering enough fuel to the SMBHs even in the earlier gas-rich 
epochs, it is hard to see how such high mass in-fluxes could be maintained in 
the local AGN. 



6. Warped accretion disc as an effective absorber 

Given unrealistically high mass outflow rates required to obscure AGN in the 
wind outflow model, I suggest that outflows, while certainly being important for 
obscuration in optically thin AGN, do not nevertheless provide the bulk of the 
obscuration that must be optically thick. This is not to say that winds are not 
important for AGN physics! 

I argue that a better physically motivated model for the absorber would 
be a warped accretion disc. The disc may be significantly warped by insta- 
bilities due t o back reaction to mass outflow, or due to the AGN radiation 
pressure (e.g. ISchandl fe Mever 19941 : iPringle 1996h . Nayakshin (2005) showed 



that an initially flat accretion disc can be quickly warped due to precession 
in a non-axisymmetric gravitational potential. The potential considered there 
was that of a stellar ring inclined with respect to the gaseous disc, a situa- 
tion which probably existed in the Galactic Centre ( Navakshin &: Cuadra 20051 : 



Paumard et al. 20061 ). Finally, there i s in general no reason for an initial disc 



configuration to be flat ( Phinnev 19891 ) 



7. Summary and Conclusions 

In this article I considered some manifestations of stellar activity in the inner 
parsecs of galaxies, inside of what is usually considered to be an AGN accretion 
disc domain. Stars are important sinks of mass during star formation and are 
important radiation and matter sources later on when stars start to burn their 
nuclear fuel. Observations of our Galactic Center highlight the problem faced by 
accretion disc models on > 0.1 pc scales. Accretion discs are very cold there, and 
one expects these discs to be making stars rather than accreting. This appears 
to be exactly what happened with two ring/disks of gas in the inner ~ 0.5 pc 
of our galaxy about 6 million years ago. If not for the star formation, Sgr A* 
could be a respectable low-luminosity AGN now (|Navakshin fc Cuadra 2003 ) . It 



seems indisputable that a working solution of the AGN/quasar feeding problem 
will have to include star formation as well. Much theoretical work remains to 
be done in this area. 
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I also considered the possible role of outflows driven away by stars popu- 
lating star forming discs like those observed in Sgr A*. It is found that these 
winds are clumpy, cold and may well contain dust, as required by the "torus" 
outflow model of Elitzur & Shlosman (2006). This model thus might work well 
for Compton-thin sources. The outflows become optically thick when they are 
powered by mass loss rates a factor of ~ 10 or more larger than the mass accre- 
tion rate corresponding to Eddington Luminosity for the AGN. This appears to 
be a general problem for obscuration by any kind of an outflow. I propose that 
warped accretion disks can provide the needed optically thick AGN obscuration. 
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